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In this paper, the issue how to introduce matter in Hofava-Lifshitz theories of gravity is ad- 
dressed. This is a key point in order to complete the proper definition of these theories and, what 
is very important, to study their possible phenomenological implications. As is well known, in 
Hofava-Lifshitz gravity the breakdown of Lorentz invariance invalidates the usual notion of mini- 
CIh, mally coupled matter. Two different approaches to bypass this problem are here described. One is 

D . based on a Kaluza-Klein reinterpretation of the 3+1 decomposition of the gravity degrees of free- 

dom, what naturally leads to a definition of a U(l) gauge symmetry and, hence, to a new type of 
minimal coupling. The other approach relies on a midi-superspace formalism and the subsequent 
. parametrization of the matter stress-energy tensor in terms of deep infrared variables. Using the 

last option, the phase space of the Hofava-Lifshitz cosmology in the presence of general matter 
couplings is studied. It is found, in particular, that the equation of state of the effective matter 
may be very different from the actual matter one, owing to the non-linear interactions which exists 
between matter and gravity. 

fL) . PACS numbers: 04.60.Bc, 04.50.Kd, 04.60.-m, 98.80.Es, 05.45.-a 

I. INTRODUCTION 

> , 

Recently, Hofava made a proposal for an ultraviolet completion of general relativity (GR), normally referred to as 
Hofava-Lifshitz (HL) gravity due to Hofava's initial inspiration on the Lifshitz theory in solid state physics. The 
salient characteristic of the HL proposal is that it seems to be renormalizable, at least at the level of power counting. 
This ultraviolet behavior is obtained by introducing irrelevant operators that explicitly break Lorentz invariance but 
ameliorate the ultraviolet divergences. On the other hand, Lorentz invariance is expected to be recovered at low 
energies, as an accidental symmetry of the theory. 

The Original HL proposal has evolved in many aspects and we count nowadays numerous sophisticated versions. 
In these, new terms have been added to the original Lagrangian, with the idea to generalize the proposal to make it 
more viable from the phenomenological perspective (see Q), and to cure the so-called strongly coupled problem (0], 
Q and references therein), via the introduction of new terms [BJ. Also, attempts to further generalize this theory at 
;_i ' the action level have been undertaken in [f| . Although, as we write this article, the consistency of the theory and 
its phenomenological implications remain still uncertain, it seems clear that the above extensions deserve a careful 
analysis. 

An important feature of the original HL theory and its modifications is the breaking of diffeomorphism invariance 
(Difn) due to the introduction of precisely those irrelevant operators that cure the UV regime. The lower number of 
symmetries in the theory, as compared to GR, produces the collateral effect that one looses the notion of "minimal 
coupling" between matter and gravity. Another equivalent point of view (at least at the classical level) comes from 
the covariant formulation of [3,0| where a Stuckelberg extra scalar degree of freedom over the metric field has to be 
introduced that may couple to matter in many different ways. In any case, for both formulations (the covariant and 
the non-covariant one) we have in principle no arguments to choose a particular type of coupling from amongst the 
most general family of couplings between the gravity and matter sectors. There is very little work in this regard in 
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the literature. For example, in Ref. Q minimal coupling is just assumed, to make more easy contact with GR. Other 
options where some particular couplings have been considered can be found, for example, in Ref. 

In the present paper, two different viable ways to approach this important problem are considered. One of them 
is a general framework that will teach us how to incorporate, in an educated manner, our ignorance on couplings 
between matter and gravity. In fact, this method can be used mostly in cosmology, but also in Black-Hole physics 
and other general situations in which certain amount of space-time symmetries are assumed. Here, we will import 
points of view and basic methodology from midi-super-space approaches and the 3+1 decomposition. The main 
idea is to parameterize the total four-dimensional energy tensor in terms of the deep IR energy variables, thus 
obtaining a formal expansion where the IR limit corresponds to the usual GR stress-energy tensor and higher-order 
terms correspond to the particular modifications introduced by the HL theory Since these IR variables represent well 
known matter that we see in our laboratories (like for example density and pressure (p,p)), they should satisfy the 
usual equation of state and conservation laws, since the theory is assumed to recover Diffi at low energies. Our other 
approach is based on a reinterpretation of the 3+1 decomposition as a form of Kaluza-Klein dimensional reduction, 
where we still have an untouched three dimensional Diffi. Then, use is made of the fact that electro-magnetic duality 
in three dimensions relates one-forms to two-forms, such that we can translate the couplings of matter with the shift 
N, into a [/(l) gauge field coupling to charged matter. At this point we recover a seance of minimal coupling based 
on the gravitational U(l) gauge theory. Obviously, this [/(l)-symmetry is only relevant for the matter sector that 
couples to N and represents, therefore, only a partial solution to the general problem. 

After defining the above frameworks, we proceed to apply these ideas to cosmological scenarios. Much research on 
this particular aspect of HL gravity has been done in the last two years [Toj]. Here, in order to address the rather 
involved issue of studying the cosmological phase space of non-minimally coupled HL gravity, we will borrow, as we 



already did in a previous work [ll| , specific techniques from the field of dynamical systems that are frequently used on 



more canonical studies applied to diverse types of cosmologies (see 1121 for a thoughtful introduction to the technique) . 
Another example of this kind of analysis of HL gravity is given in [13j . Later in this paper, we specifically study the 
cosmological phase space of the HL model with our new matter couplings. Our current investigation is focused on the 
introduction of non-standard couplin gs b etween matter and gravity and should be understood as a genuine extension 
of our previous work on this subject [ll| to cases of real phenomenological interest. 

The paper is organized as follows. In Sect. |TT] we present a short overview of the relevant modification of the HL 
theory we are considering, namely the inclusion of a minimal potential defined in 0], and the so called "healthy 
extension" of Q, that presents a whole family of new terms related to the lapse function N. In Sect. Mil we present 
our framework to study matter couplings while in Sect. IIVI wc describe an application of the above frameworks, to 
cosmological scenarios with our generalized matter couplings. We are then able to characterize the classical phase 
space, discussing its structure in depth, in particular all the fixed points and their nature as repellers or attractors in 
the theory with matter couplings, in each of the corresponding cases. Finally, in Sect. [V] we summarize the results 
obtained, giving some perspectives for further work. 



II. HORAVA-LIFSHITZ WITHOUT MATTER 



In HL gravity, the gravitational dynamical variables arc defined to be the lapse N, the shift iVj and the space metric 
gij, Latin indices running from 1 to 3. The space-time metric is defined using the ADM slicing of space-time, as 

ds 2 = h^dx^dx" = -N 2 dt 2 + gij (dx l + N i )(dx j + N j ) , (1) 

where N l = g ZJ Nj, as usual. The action S is written in terms of geometric objects, covariant under 3d-diffcomorphisms, 
characteristics of the ADM construction, like the 3d-covariant derivative Vj, the spatial curvature tensor R%jki, and 
the extrinsic curvature Kij. They are defined as follows, 

Wjki = W*ji,k - W*jk,i + W m 3 iW\ m - W n jk W i im , (2) 
where W l ji are the Christoffel symbols (symmetric in the lower indices), given by 

w ji = 7;9 im (9jm,i + g m i,j - 9ji, m ) ■ (3) 
The Ricci tensor is obtained by contracting the first and the third indices 

Rij = 9 kl R lkj i and R = gVRij , (4) 
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while the extrinsic curvature is defined as 

Kij = ^(~9ij + Vity + VjNi) , (5) 

where the dot stands for time derivative. 

In terms of the above tensor fields, the HL action can be written as 

S = J dt dx 3 Ny/g (C klneUc - C 

potential H~ ^matter ■) j (6) 

being the kinetic term universally given by 

Ckinetic = "fA , ,A" - \K 2 ) , (7) 

and with a and A playing the role of coupling constants. Originally, the potential term was a generic function of Rijki 
and Vj but in Ref. p it was realized that this generic function should also depend on a.j = V,; ln(iV). 

The action generically breaks covariance down to the subgroup of 3-dimensional diffeomorphisms and time 
reparametrization, i.e. x — > x(t,x) and t — > t(t). Assigning dimension -1 to space and dimension -3 to time, it 
can be seen that it is enough to restrict the potential to be made out of operators up to dimension 6, in order to get 
a power-counting renormalizable theory. 

Here, we will work with a potential which corresponds to the more general choice available, composed by the 
potential defined in [H (the SVW case), which depends on Rijki and Vj, 

^potential — SVW P&iRjHVB?" + (3 7 RV 2 R + ijriiilK + f3 5 R(R 3k W k ) (8) 
+/3 4 i? 3 + (3 3 R ]k R ]k + (3 2 R 2 + PiR + (3 , 

and with the addition of all the general terms as suggested in [f|. We collect all these terms in an implicit form 

£potential-a,i = ^ Tn O n (di, Vj, Rijki), (9) 

where O n are general operators of maximum dimension 6 and 7„ the corresponding coupling constants. 

In our calculation, we have worked out all the independent terms of these operators and have chosen to display 
only, as the representative operator for each class, the ones with less derivatives acting on a single at 1 . We therefore 
obtain: 

At order 2 

70^ + 71 a 2 , (10) 



at order 4 
and at order 6 



73 a 4 + 74 a 2 D % a % + 75 D i a i D j a j + 7 6 a 2 R + 77 r/V//, , + 7 8 D 2 R + 79 //,,/," ' + 710 R 2 , (11) 



7x1 a 6 + 7 12 a 4 R + j 13 ^a.ajR^ + 7 14 a 2 R 2 + -y 15 a.a^R^R + 7i 6 a t a k R 13 R) + 717 R 3 

+7i8 R kl R l3 Rj + 719 RR ij Ri 3 + 720 "-/."''/.',, + 721 a A D t a l + 722 a 2 D ia l B 3 a? + 723 D i a i D j a' ' D k a k 

+724 />,</' /V/'/W + 725 D. l a j D k a t D k a j + 725 a 2 a k D k R + 727 a t a k D k D t R + 7 28 a 2 D 2 R 

+729 D 2 D 2 R + 730 a i alR tj D k a h + 731 a t DR l D k a k + 732 a'R tk D t R k + 733 IV li ,,,!)■ I!) 

+734 RD 2 R + 735 D*RDtR + -3,, /."'/ >,/>.,/.' . (12) 

The above list includes previous terms of the SW-potential and a number of new terms due to the appearance of the 
new field a. Notice that we have defined the potentials such that they have to be multiplied by N^Jg, which explicitly 
contains N. Therefore, even though some of the new terms in the potential do not exhibit an explicit coupling with 
aj, they are not equivalent to any term of the SVW-potential. 

At this point, other phenomenological and theoretical considerations may help us constraint the range of values 
the different couplings should take. For example, in Ref. [lH it was found that ghost instabilities are present if 



1 This is an arbitrary basis that, nevertheless, fixes our conventions. 
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A e (1/3, 1), that the cosmological constant is negative for the detailed balance potential, a > 0, and so on. Here, we 
will constraint as little as possible the different ranges of values on each coupling constant to see how much information 
comes out of the dynamical system approach itself. Then, we will add this information to the constraints arising from 
other considerations, to finally obtain the most promising form of the potential. In particular, we will take A different 
from 1/3 (corresponding to the scale invariant case) as the only limitation on its range. 

III. MATER COUPLINGS 

Once we have explicitly defined the extension of the HL theory we will be working with — at least what concerns its 
gravity sector — it is time to focus now on how matter is to be coupled to gravity. As we mention in the introduction, 
due to the reduction of the symmetries present in the theory, we have no longer a valid argument to define a minimal 
coupling. In fact, this is more dangerous than what it may naively seem, since for example, different particles will 
have in general different dispersion relations, depending on their couplings with the gravity sector; also, since we have 
more geometric invariants, there are many more ways to construct couplings to matter. We postpone this line of 
thought to future work, to focuss here just on the the simplest problem of parameterizing the possible form of all 
these different couplings in terms of physical quantities. 

A. Midi-superspace approach 

The first approach is based on two main assumptions. First, that in the deep IR regime one should recover 
diffcomorphism invariance and hence, that the corresponding IR stress energy tensor T Ml/ should be divergence-free. 
Using the above, we can formally expand the full stress energy tensor, 7j w , in terms of physical observables, defined 
in the deep IR regime alone. Second, we assume that we have at disposal a considerable number of symmetries, like 
in cosmological models or Black Hole physics, which allow us to write 7^ v in terms of just a few variables, as the 
energy density p, the pressure p, fluid velocity u, etc. The above assumptions imply that we can write 

% v = T^(p,p, ...)+ AT^(p,p, . . .), (13) 

where AT /JJ/ (p,p, . . .) is the leftover contribution, made out of irrelevant operators and other such terms which will 
anyway decouple at low energies. The above expansion can be understood as a sort of derivative expansion in the 
gravitational coupling with matter fields. From the point of view of the covariant formulation, where a Stuckelberg 
field 4> is added on top of the metric g, what we are here doing is to separate the matter Lagrangian C ma tteri^\ where 
ip represents matter fields, into a part which is minimally coupled L rn i n {g;il)) and the rest of it, £ non - m i n ((j), g;if)), 
which generically has couplings to and g, e.g., 

^matter 

(4') = C min (^;g) + C «>,<?; 0)- (14) 

Using the above ideas, we can describe our HL theory in terms of a Lagrangian formulation where the gravity 
sector is given in terms of midi-superspace variables, while the matter sector is described in terms of hydrodynamic 
variables, like (p,p). For example, consider the case of the celebrated FRW ansatz. Here, due to the symmetries 
imposed, 7^ depends on two variables (p,p) only and it can be written so that it is diagonal. Then, our previous 
considerations translate into the following expansion 

T oa =p T = p + S Q p + S lP 2 + S 2 Rp + S 3 Rp 2 + S 4 R 2 p + S 5 p 3 + S 6 p, (15) 

— jj-^- = P T = p + VoP + mp 2 + mRp + vaRp 2 + V±R 2 P + S 5 p 3 + VeP, (16) 

where, for consistency with the gravity sector, we have limited the expansion to operators up to order 6. Note that 
here (p,p) represent standard matter and, therefore, they fulfill the usual linear relations 

p = wp . (17) 

Moreover, owing to the non-minimal couplings to gravity, the above equation of state gives rise to non-linear relations 
for the effective total energy density and pressure (p T ,p T ). These equations, together with the equation of state, will 
define the type of fluid we can consider in our HL cosmology where matter is coupled to gravity in the most generic 
form. 
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Coming back to the general case, as in GR, minimization of the action S, upon variation of the metric in the pure 
gravity sector, defines the two-index tensor HL UU1 

HL^v = ^ — j ^ — dtdx 3 Ny/g (£ kinetic ~ ^potential)^ ■ (18) 

This tensor can also be decomposed into an IR part, corresponding precisely to the Einstein tensor G U u, plus a 
leftover, characteristic of the HL theory, which we write as 

a F = ^ (G M „ + AHpu) , (19) 

where k 2 is the gravitational coupling constant (in natural units k 2 — 8itGn)- Therefore, the form of the field 
equations obtained by minimizing the HL action with respect to the metric is 

Guv + A% = k 2 {T nv + AT nu ) . (20) 

At this point we still have to add a last constraint, that comes from taking the 4d divergence to the gravitational field 
equations (|20p and using the usual Bianchi identities for G u „ and v . This yields namely 

4V ^ (AHLftv) = 4V (ATf,, u ) , (21) 

where 4V is the 4d covariant derivative. It is important to recall that these equations have to be satisfied only 
on-shell, since they actually come from the field equations. 

Summarizing, after the whole derivation has been carried out, our final set of equations, that define our recipe to 
deal with general couplings of matter to HL gravity, is given by Eqs. (|17I20I21[) . Observe that the above expressions 
are written in terms of a set of operators of order less or equal than 6, and are made out of geometric objects, like the 
3d curvature R, the vector aj, and the 3d covariant derivative V, and of IR hydrodynamic variables, as for example 
(p, p)- 



B. U(l) gravitational coupling 

As we have argued previously, due to the break down of diffcomorphism invariance, there is no more a clean 
argument to define what it will be a minimal coupling between gravity and the matter sectors. One could still try, of 
course, to use a principle based merely on simplicity, but even then, it is a fact that such a principle will by no means 
be universal, not to talk on its grounding from pure physical considerations. Instead of following the above line of 
thought, we will here argue that the remaining symmetries of the theory are still strong enough in order to help us 
find a clear guiding principle which, under some general assumptions, will deliver a well defined definition of what a 
minimal coupling in HL theories of gravity should be. 

To illustrate our idea, let us consider the more relevant part of the potential terms in the gravity sector of the HL 
Lagrangian, 

J dtdx 3 Ny/g (a R + a ia l a z ) . (22) 

We can always device a conformal transformation in the 3d metric, of the form g^ — N~ 2 gij, such that, in the new 
frame, the vector a and the 3d curvature tensor R arc canonically normalized with respect to the 3d metric (jij, i.e. 

dtdx 3 yjg~ (a R + -//"(/,« , ) , (23) 

where 7 = (ai — 2ao) and where we have discarded pure boundary terms for simplicity. In this new frame, is a 
3d one-form, which can be transformed into a two-form Fjk via Hodge-duality. The resulting two-form is naturally 
described in terms of a one-form gauge potential bi, what immediately leads to a gravitational U(l) gauge symmetry: 

" e ijk F jk , Eij = dibj - djbi, (24) 



2 3 

where e is the 3d Levi-Civita pseudo-tensor in the new frame. The constant factor in the above definition is chosen 
to have a canonically normalized kinetic term in the corresponding action, 

dtdx 3 y/l (a R + -l^g^g kl F lk F^ , (25) 
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where g~ 2 = 7. 

From this point of view, we have a local U(l) symmetry that may be used to define the meaning of our new 
"minimal coupling" with this gravity sector. Notice that, in order to respect such symmetry, we can only write gauge 
invariant operators in our matter Lagrangian and, therefore, it is clear that our minimal coupling should be written 
in terms of covariant U(l) derivatives, of the form 

A = V, - ib t . (26) 
For example, the Lagrangian for a charged scalar field, tp, minimally coupled to the gravity sector, will be given by 

1 



£ = (Vs) 



(27) 



where, as a complementary principle, we have used the fact that, at low energy, we should necessarily recover the 
usual meaning of minimal coupling. 

The case of vector fields and their possible couplings to gravity is more subtle. A gauge field A can be coupled to 
a conserved current J in the minimal form as J 1 A,;. In 3d we have a topologically conserved charge, made out of our 
gravitational gauge field b, namely J 1 = e i3 Fj}.. This current is conserved by the corresponding 3d Bianchi identities 
and, therefore, opens the possibility of including such a term in the usual gauge invariant action, producing thereby 
a sort of Chcrn-Simons term coupling the matter gauge field to the gravity gauge field. 

The above idea is far from having been investigated in detail and deserves a lot more attention. We believe that it 
is indeed an interesting starting point to study such a complicated subject as is the coupling of matter to HL gravity. 
Note that, as an outcome of these ideas, we are left with a sort of mixed argument, based on symmetries and the low 
energy emergence of 4d diffeomorphism invariance, which still does not have a unique, well defined, general meaning 
as a "minimal coupling" in HL gravity. Nevertheless, we have clearly achieved some relevant improvement, since at 
least the gravity sector related to all the operators made out of a.; is under control, due to the gravitational U{1) 
symmetry. It will be very interesting to study the implications of the above ideas to BH physics, where the form of the 
metric clearly induces the coupling with the vector. Hopefully these studies will appear soon elsewhere. In what 
follows here, we will directly concentrate on cosmological scenarios where, owing to the form of the typical ansatz, 
the C/(l) gravitational symmetry is trivial and we still do not take real advantage of the above mechanism. 



IV. NON-MINIMAL COSMOLOGY 



In this section we apply the above methods to study the impact of non-minimal couplings 2 to HL gravity in 
cosmological scenarios. We start, as our basic ansatz, from the usual FRW metric corresponding to homogeneous and 
isotropic space-time, i.e. 

ds 2 = -dt 2 + aitf^dx'dx 1 , (28) 
where 7^ is a maximally symmetric metric of constant curvature k = (—1,0, 1), and a is the scale factor such that 

Rij = —^gij j Kij = ~Hgij , H = — . (29) 
With these symmetries, the total stress energy tensor T^v can be written as 

T = p T dt ®dt + p T g l3 dx l ® dx> , (30) 
where p T and p T are the effective total energy density and effective total pressure, which can be expanded as in 

Eq. cas). 

The gravitational sector is drastically simplified in the FRW ansatz where only the terms corresponding to the 
coupling constant /?o, Pi, P2, P3, Pi, P5, and Pe in Eq. ([5]) are not identically zero. In particular, this means that all 
terms related to the field d\n(N) do not contribute. Following, the work of Q, we assume that our coordinates are 
such that c = 1. Wc have also used time re-parametrization plus three-dimensional diffeomorphism invariance, in 
order to eliminate redundant degrees of freedom. 



2 We use this term, in contraposition to the usual minimal coupling in GR, to denote a general coupling of matter to HL gravity. 
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In this setting, the non-zero components of the HL tensor are 



HLij ■ 



3 [ 



2 



A- 
2H + 3H 2 



X3k 2 
2a 4 
k 



2a 6 



a z 



,Y3 



fc 2 



6a 4 



X4k \ 
2a 6 J 



(31) 
(32) 



where we set 



16nG N = i , A 

X3 = 12a(3/3 2 +/3 3 )a 2 , 



, 1 = -Pi* 2 , 
Xi = 24 (9ft + 3/3 5 + /3 6 ) a 4 



Therefore, we have 



AHL 00 

Aft,, : 



X3k 2 



2a 4 

£ ( 2ff + 3ff 2 



2a 6 J ' 
6a 4 



2a 6 



while the constraint equation (|21[) in this framework can be written as 

1 



(AiZLoo) + 3ff f Ai^Loo + | AJ 



(AToo) + 3# AToo + ^AT- 



(33) 

(34) 
(35) 

(36) 



It is possible to understand this equation as a measure of up to which level the failure of diffcomorphism invariance 
in the gravity sector is transmitted to the matter sector. Nevertheless, for the FRW ansatz, it is not difficult to verify 
that 4V (AHL^v) = 0, due to the homogeneity of the ansatz and, therefore, we are left with the following pure 
constraint on the matter sector: 



AToo + \&T 



0. 



(37) 



This non-trivial relation for the matter sector ultimately reduces the number of independent couplings that characterize 
the expansion (|15H16[) . In fact, one finds the following relation between these couplings: 



S 6 + 5 W = ?/o + T)6W , 8i(l + 2w) = 7]i , 

8 5 {2 + 3w) = 775, 8 2 {w + 2/3) = m , 

<S 3 (5/3 + 2 w ) =773, S 4 (w + 4/3) = 774, 

for general values of w. 

The field equations for the gravity sector, HL^ V = \T^ V result into the following expression 

1 - 2i\ H 2 - Mi - XL - Xik 2 _ Xik _ 2 n T _ n 

1 2 / n 6a 2 6 6a 4 6a 6 K P ~ U ' 



(38) 



(39) 



which, together with Eqs. (|17j) . (|15H16j) . and (|38|l define our cosmological system. The presence of non minimal 
couplings of the matter terms conveys the idea that, because of the Lorentz violation, what drives here the cosmological 
expansion is no more p, but actually p T . As a consequence, depending on the values of the parameters 8, a certain type 
of matter can rather behave effectively as a fluid with different thermodynamical properties. The natural question 
being then, wether such behavior can possibly help explaining some important features of the observed Universe. In 
order to answer it, let us first investigate Eq. (|39[) for the two basic classical matter types: dust and radiation. 



A. The dust case (w = 0) 

In the case of dust, the last equation above reduces to 



(40) 
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which is exactly the same system treated in with the only difference that the coupling constant is now modified 
by the parameter Sq. That is to say, we now have control on the way in which matter couples to gravity at the 
cosmological level. Note however that such small change can induce a great deal of difference in the behavior of the 
respective cosmology. For example, in the case 5$ < — I, p T has a negative coupling constant, which generates an 
effective dark energy, even if the standard matter has the usual thermodynamical properties. In standard GR this 
would lead to an irreparable inconsistency of the theory. However, in HL gravity the presence of the additional terms 
does not exclude such case. Of course, once the coupling constant is set to be negative it stays so forever, and its 
effect on the cosmic processes typical of the dust era should be investigated carefully. 



B. The radiation case (w = ~) 

Let us consider now the case of radiation. The cosmological equations read 



tt2 I i 3A 5 5 p 3 2kS 3 p0 2 ft ft ft k X2 Xi n (AU 



3a 10 


9a s 


6a 6 


6a 4 


6a 2 


6 


H 2 - 




8k5 3 p 2 


h 


P2 


As 


27a 12 


3a 10 


3a 8 


a 6 


6a 4 



2 J V 2 J 27a 12 3a 10 3a 8 a 6 6a 4 6 



being 



ft=p (l2 2 ^4 + <5 1(0o ), (43) 
h = k 3 X i + I2k6 2 p , (44) 
ft = k 2 X3 + 2(5 Q + 35 6 + 3)p . (45) 

Differently from the previous situation, here the non-vanishing pressure "switches on" new terms associated to the 
additional matter couplings. Looking at the general structure of the equations above, however, it is clear that the 
effects of these couplings will only be relevant at early times, namely when the scale factor a is particularly small. 
Thus, we can conclude that the introduction of a full matter coupling will influence the evolution of the early Universe, 
mainly. 

It is also interesting to note that, when radiation dominates, the presence of this coupling induces differences 
between the HL cosmology and the GR one, even in the case of spatially flat solutions. Such difference is not present 
in the dust case. 

In order to have a more clear idea of the effects of the matter coupling in Eqs. (|41j|42p we can use the dynamical 
system approach. Given the high number of degrees of freedom of the system, we will limit ourselves to consider the 
finite analysis only. Following the method of [ll| , we define the variables 

X = Xl Y = 265P " Z = 4k63p2 ° 



3(3£-2)iT 2 ' 27(3£ - 2)H 2 a 12 ' 3(3£ - 2)H 2 a 10 ' 

D Pi q P2 rp Ps ,, R \ 



K 



9(3£ - 2)H 2 a 8 ' 3(3C-2) J ff 2 a 6 ' 3(3£ - 2)# 2 a 4 ' 

3(3£ - 2)ff 2 a 2 : 



through which we will characterize the phase space, and a logarithmic time N = In a. 
Then, the resulting dynamical system is 
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-25 4 
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-5Y- 


4 4Z4 




o = 


= 1 + K + R 


+ 54 


-T + 


n 


-Y + 







(47) 



(48) 
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where the "prime" indicates derivative with respect to N. The structure of the system reveals that the phase space 
is divided into different sectors, delimited by invariant submanifolds, and that, as a consequence, no global attractor 
can actually exist 3 . Therefore, any orbit that can have physical interest will be realized by only using a restricted set 
of initial conditions. 

The fixed points can be found, as usual, by setting the lhs of the equations equal to zero. Then, the solutions 
associated to the fixed points can be found using the general expressions 

H = aH 2 , (49) 

P „ 3 , 

- =4H = -, (50 

P a{t-t )' v ' 

a = -1 + 3R t + 2S t + Ti - X, + 5Y t + 4Z t , 

where the "i" subscript represents the value of the corresponding variable at the fixed point. Subsequently, using the 
Hartman-Grofimann theorem, we can investigate their stability. The finite fixed points, their associated solution and 
their corresponding stability are all summarized in Table [I] 

As expected, we find fixed points associated to the dominance of the different terms in Eqs. (|4HI42[) . These points 
are characterized by the corresponding expansion rates of the scale factor and of the dissipation of the energy density. 
Choosing initial conditions in which all the variables are negative (which correspond to a specific set of constraint for 
the coupling constants \i an d 6i ) , we will select the sector of the phase space in which the orbits have access to all 
fixed points. In the same way, cosmic histories that surely avoid any of the states described by one of such points can 
be selected by modifying the values of the couplings. 

Note also that none of the fixed points related to the matter couplings appears to be characterized by an accelerated 
expansion, regardless of the relative sign of the coefficients of those terms. We will discuss in the conclusions the 
possible consequences of the presence and nature of these points. And note moreover that, in spite of all these changes 
one still has a de Sitter attractor that is generated by the cosmological terms in the gravitational part of the action. 
Therefore, the introduction of the matter couplings does not, in principle, compromise the formation of a dark energy 
era, in spite of the nature of the cosmic fluid and of the structure of the non-minimal couplings. This is a quite 
remarkable result. 

Substitution in the cosmological equations reveals that the fixed points A, C — Q are not actual solutions of the 
cosmological equations. However, this is not a serious problem, since all these points are unstable and, therefore, they 
may only correspond to an approximation of the general solution represented by the full orbits. On the contrary, the 
point B corresponds, instead, to an exact solution of the theory. 



To be precise, this role could be taken by a point at the origin. However, as we will see, such point would be always unstable. 
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TABLE I: Coordinates, solutions, and stability of all the fixed points of the dynamical system (1471) . 



Point Coordinates Solution Energy Density k Stability 
[X,Y,Z,R,S,T] 
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saddle 
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saddle 


Q 


[0,0,0,0,0,-1] 


a 


= ao(t — 


to) 1/2 


p 


= it- 


to)' 


-2 





saddle 



V. DISCUSSION AND CONCLUSIONS 

We have proposed in this paper two different ways to study the coupling of HL theories of gravity with matter. 
First, we have devised a very natural procedure to study the effect of matter couplings in a Hofava-Lifshitz theory of 
gravity, based on the imposition of the Bianchi identities on a midi-superspace approach together with the assumption 
of a well-behaved IR behavior. The use of these geometric relations has revealed that, when matter is present, we 
need to supplement the theory with an additional constraint that ensures the consistency of the gravitational and 
matter sectors. Second, we have introduced an alternative definition of 'minimal coupling' based on the U(l) gauge 
symmetry present in a very important sector of the theory. In this way, we are able to couple the matter fields via a 
gauge invariant formalism. These couplings are both relevant for BH theories and also for perturbation theory in a 
cosmological scenario, but irrelevant for the study of cosmological solutions, since in this last case the corresponding 
gravity sector is identically zero. 

Regarding the application of the above ideas to cosmological scenarios, in order to put a first set of constraints on this 
new version of the theory we have chosen to analyze its Friedmannian cosmology. The high symmetry of cosmological 
space-times greatly simplifies the equations and allowed us to introduce, in a simple way, a set of (arbitrary) additional 
matter couplings. In particular, the constraint simplifies into algebraic relations among the coupling constants of the 
additional matter terms, which can be used to understand their specific physical role. 

A first interesting finding here was that, because of the nature of the constraint, the behavior of the HL Universe is 
very sensitive to the thermodynamical properties of matter. The cosmological equations arc enriched with additional 
terms and both their structure and their number change with the barotropic factor. This implies that the resulting 
phenomenology is expected to be quite different from one case to another, when one considers Universes dominated 
by different types of matter. 

In addition, although in principle the new matter couplings influence the entire cosmic history, from the structure 
of the equations it was easy to conclude that the new terms have more weight at early times. This suggests that 
the most important deviations from standard GR will be evidenced at an early epoch only. This has an important 
consequence in what concerns the testability of the theory, because the early history of the Universe is very tightly 
constrained by several independent sources of evidence. 

In order to examine more in detail the specific phenomenology, we have used dynamical system techniques. In the 
case of dust we have found that the resulting model is basically the one already analyzed in [ll| , with the difference 
that now the sign and the strength of the matter couplings become free parameters or the theory. However, in the 
end this fact does not have any influences on the number of the fixed points nor on the solutions associated to them, 
but only on the set of initial conditions that lead to a certain behavior. 
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The case of radiation is more complex because of the additional terms that appear in the cosmological equations. 
In phase space such terms are associated to additional fixed points (.4, C — Q). If all variables are negative, one can 
obtain the richest possible behavior for the corresponding cosmology, i.e., orbits that are able to "touch" the entire set 
of fixed points (albeit they are all unstable). The presence of such fixed points and, in particular, the existence of the 
solutions associated to them, constitutes one of the main results of the present paper. In fact these fixed points allow 
us to draw quantitative conclusions on some of the best studied phenomena of the early Universe, as nucleosynthesis 
or cosmic microwave background (CMB) physics. 

It is natural to ask if one can choose the above mentioned constants is such a way that the models develop and 
early-time accelerated expansion phase which might work as inflation. The answer to that is, unfortunately, negative, 
for the time being, and the associated dynamical system gives a clear way to see it. In fact, one can prove that, given 
a generalized Friedman equation of the form 

^ «L + _5*. + ... + .EL (51) 

a mi a m 2 a mi ^ > 

only terms in which < mj < 2 can generate fixed points which may represent accelerated expansion. This is clearly 
not the case for any of the terms appearing in (|40|) or in (|41[) . What means that, in the phase space, there is no fixed 
point which could correspond to an inflationary era. However, the situation is not as bad as it might seem at first 
sight. It has been already suggested that, within the HL framework, other mechanisms might substitute an inflationary 
phase @. Our results seem to support this last, alternative scenario instead of the classical one. Therefore, in spite 
of the degrees of freedom added by the additional matter couplings it is very difficult to reproduce any accelerated 
expansion and the only way in which one can achieve a dark energy era is to introduce the cosmological constant "by 
hand" in the HL potential. This appears as a clear limitation of the theory since it does not seems to provided a 
natural solution to the quantum gravity phenomenon of dark energy 4 . 

Anyway, the important issue to be remarked is that, thanks to a consistent introduction of the matter field, the HL 
theory in this new form finally becomes directly testable against some of the most accurate data we posses at present. 
And these data hold the promise to be able to verify directly its validity. 
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